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Abstract 

We compute the matching conditions for a general thick codimension 2 brane, 
a necessary previous step towards the investigation of gravitational phenomena 
in codimension 2 braneworlds. We show that, provided the brane is weakly 
curved, they are specified by the integral in the extra dimensions of the brane 
energy-momentum, independently of its detailed internal structure. These gen- 
eral matching conditions can then be used as boundary conditions for the bulk 
solution. By evaluating Einstein equations at the brane boundary we are able 
to write an evolution equation for the induced metric on the brane depending 
only on physical brane parameters and the bulk energy-momentum tensor. We 
particularise to a cosmological metric and show that a realistic cosmology can 
be obtained in the simplest case of having just a non-zero cosmological con- 
stant in the bulk. We point out several parallelisms between this case and the 
codimension 1 brane worlds in an AdS space. 



*E-mail: ignacio .navarroQdurhain. ac .uk 
^^E-mail: jsantiag@fnal.gov 



1 Introduction 



New cosmological observations seem to imply that the expansion of our universe is 
currently accelerating [1] , driven by a dominant component of the energy-momentum 
tensor (EMT) with an equation-of-state parameter w close to -1 (the so called dark 
energy). The observations have made the cosmological constant problem a very pressing 
one: to the traditional puzzle of an (almost) vanishing vacuum energy now cosmologists 
(and particle physicists) wonder why is its magnitude comparable to the matter energy 
density today. Recent analysis of the data [2] point to an even more bizarre situation: 
the best fit to observations agrees with a dark energy equation of state with a; < — 1. So 
the family of problems associated with the vacuum energy (the cosmological constant 
problem and its smaller cousin, the coincidence problem) could grow in the near future 
with a new member: why is the vacuum energy growing in time? The problem with 
this possibility is that it is not easily accommodated in generally covariant theories as 
long as the matter EMT satisfies the usual energy conditions [3] . In the same fashion, 
other existing observations also suggest modified gravitational dynamics^. 

Weinberg's theorem [6] shows that standard approaches (by which we mean 4D field 
theories based on General Relativity) to the cosmological constant problem are very 
likely to fail and therefore more exotic ones should be tried. In particular having more 
than four dimensions in a Kaluza-Klein fashion does not seem to improve the situation, 
since the extra dimensions are small and gravity is effectively four dimensional below 
some scale. In this effective theory one will face the same problems as in any 4D theory. 
Thus, within this class of theories, anthropic considerations seem at the moment the 
only framework capable of explaining some of the large scale properties of our universe 
[7]. Brane- world gravity, on the other hand, does not belong to this class of theories, 
since it is not guaranteed that the low energy description of gravity can be obtained 
from a generally covariant four- dimensional Lagrangian, i.e., there is not necessarily a 
four dimensional description of the gravitational sector. In these models one assumes 
that the Standard Model fields are confined to some submanifold of the whole space- 
time. One can think on the Standard Model fields as the zero modes of topological 
defects of higher dimensional field theories [8] or the gauge theories hving on the world- 

^One can mention the flatness of galaxy rotation curves, that can be explained using the dark 

matter hypothesis but that can also be regarded as pointing towards modified gravitational dynamics 
[4], or the more nearby measured Pioneer 10/11 anomalies [5]. 
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volume of the string theory D-branes [9] . Fermionic fields and gauge interactions can be 
in this way clearly lower dimensional but to elucidate if the gravitational interactions 
can be well approximated by 4D Einstein gravity is not so obvious. The hope is to 
find in this context a theory that shares the good conceptual advantages of 4D General 
Relativity (gravity as geometry, background independence...) but can yield a realistic 
but non-standard cosmology (or gravitational dynamics). 

The cosmology of codimension one braneworlds is quite well understood. It is pos- 
sible to recover something close to standard cosmology with corrections that take the 
form of "dark radiation" plus terms involving the matter energy-momentum tensor 
squared [10,11]. Although the situation is not as good in the understanding of codi- 
mension two brane worlds, great progress has been made recently in their investigation. 
The fact that one can find solutions for a fiat brane in a given setup for any value of 
its tension has encouraged many authors to try and attack the cosmological constant 
problem in such scenario [12-17]. The effect of the brane tension in these models is 
simply to produce a deficit angle in the transverse space, without further implications 
for its induced metric. The two dimensional space transverse to the brane acquires lo- 
cally the geometry of a cone, with the brane situated on its tip. However one problem 
of codimension two brane-worlds is that with a deficit angle one can only generate a 
two dimensional delta function in the Einstein tensor that is proportional to the brane 
induced metric. This means that one can only find nonsingular solutions if the brane 
EMT is proportional to its induced metric, i.e. it is pure tension [14,15]. Thus the 
solutions found in [12] cannot be extended to a general brane EMT in the thin brane 
limit. This limit is indeed singular for a general brane and, as such, makes all the ar- 
guments about the nature of gravity (and self-tuning) in codimension two braneworlds 
in Einstein gravity dodgy when working with 5-like branes. To make things worse, the 
very Einstein equations imply that in the case of an infinitesimally thin pure tension 
brane, the deficit angle is space and time independent. This situation is very similar 
to the cosmic string models studied in 4D (see [18] for a rigorous treatment of codi- 
mension one and two sources in 4D General Relativity). It is therefore not sensible to 
ask questions like what happens if there is a sudden phase transition that changes the 
tension of the brane as such thing is not allowed by the equations of motion in the thin 
brane limit. In other words, self-tuning has to be formulated as a dynamical process 
and delta-like codimension two branes do not allow to do that. In [15] a possible way 



2 



out of this situation was proposed by adding the Gauss-Bonnet term to the Einstein- 
Hilbert Lagrangian. In this case, the thin brane hmit is well defined and one finds the 
remarkable result that four dimensional Einstein gravity is recovered as the dynamics 
for the induced metric on the infinitesimally thin brane. (Similar conclusions have 
been obtained in [19] at the linearised level.) See also [20,21] for a different approach 
to codimension 2 braneworlds in Gauss-Bonnet gravity. 

Motivated by the previous considerations, we abandon in this paper the thin brane 
idealization, and compute the "matching" conditions for a general thick codimension 
two brane in Einstein gravity^. By matching conditions we mean equations that relate 
the values of the first derivatives of the metric (with respect to the orthogonal coordi- 
nate, r) at the brane boundary with the brane EMT. In fact, since we are dealing with 
a thick defect (with a singular thin limit), it is not clear that one should be able to do 
this without knowing all the precise small scale structure of the brane. If we need this 
information in order to find the solution, the matching condition approach would be 
of no use, since one should solve the equations independently for different microscopic 
brane models, and one could obtain different results for different models even if the 
total energy-momentum carried by the brane is the same. Gravitational physics would 
then in general depend on the ultraviolet details of our theory and therefore model 
independent assertions would be difficult to make. We will see, however, that one can 
obtain this set of equations depending only on the integral of the brane EMT along 
the extra dimensions as long as the parallel derivatives {i.e. with respect to the brane 
coordinates) of the metric, and in particular the Ricci tensor of the brane induced 
metric, are small enough: in this case our matching conditions do not depend on the 
inner structure of the brane. How small is "small enough" will be made clear in the 
next section but one can argue that this situation is quite general in the sense that it 
is natural that the presence of the brane induces much larger gradients in the radial 
(transverse) direction that in the longitudinal ones. In particular this is clearly the 
case if we are interested in late time cosmology. Of course that does not mean that 
nothing can be said about very early cosmological times or other situations with strong 
gravity effects, but one should keep in mind when dealing with such situations that 

^We do not introduce the Gauss-Bonnet term now because once one considers a thick brane the 
Gauss-Bonnet contributions, although crucial in the thin limit for obtaining a regular geometry, will 
be subleading unless the Gauss-Bonnet coupling is very large or the brane is extremely thin. 
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our matching conditions imply certain approximations that break down when the 4D 
curvature is very large. Once we consider a thick brane, departure from pure tension 
is allowed and the deficit angle can develop a space-time dependence, thus questions 
about the self-tuning behaviour of the system are again legitimate. 

In the next section we will carefully explain our assumptions and approximations 
and we will obtain the equations that relate the total brane energy-momentum with 
the deformation of spacetime it causes, the so called matching conditions. We will 
specialize these equations to the cosmological case in the third section. Evaluating 
Einstein equations just outside the brane and using the matching conditions we will be 
able to obtain the equations that govern the cosmological evolution of this braneworld. 
In this respect, we follow a procedure completely analogous to that of the codimension 
one case [11]. 

2 Matching conditions for a codimension 2 brane 

In this section we will try to answer the following questions: what is the effect on space- 
time of an energy-momentum distribution that can be interpreted as a codimension 2 
defect in six-dimensional Einstein gravity? Is there a way to characterise this effect 
without knowing all the precise small-scale structure of the brane? We will see that 
the answer to the second question is yes, provided certain conditions are met, and we 
will also provide a (partial) answer to the first question. The needed conditions have 
the interpretation of requiring a weakly curved brane. 

In this section we will consider the following quite general ansatz for the metric. 



where, as usual, denote four non-compact dimensions (including a time-hke one), 
II — 0, ... ,3, whereas r, 9 denote the radial and angular coordinates of the (compact 
or not) two extra dimensions. This means that in particular the following boundary 
conditions hold (in order to avoid singularities) at r = 



where a prime denotes derivative with respect to r. We have assumed a rotational 
symmetry around the codimension 2 submanifold defined by the condition r — 0. The 



ds'^ = g^^^ix, r)dx''dx'' - dr'^ - L{x, rfd9^, 



(1) 



L{x, 0) = 0, L'{x, 0) = 1, drg^^ix, 0) = 0, 



(2) 
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metric is determined by the Einstein equation, that can be written as 

M^R^J = T^- (3) 

where is the 6D fundamental mass, Tmn is the EMT and T = T^^ its trace. 
The brane will be a cylindrically symmetric extended object that fills the region with 
r < e. Since we are trying to deal here with a general situation, and we do not have a 
particular microscopic theory for the brane, we cannot provide a precise definition of 
the brane width parameter, e, but in a given particular model it should not be too hard 
to provide a strict definition of it. In any case, the same results should be obtained 
taking a different definition of the brane width {i.e a different splitting into brane- 
bulk of the total 6D EMT), as long as we are considering the same energy-momentum 
distribution. 

The presence of the brane induces a strong r-dependence of the curvature tensor 
that should be reflected on large r-derivatives of the metric. It is natural then to 
assume that the brane has the effect of producing mainly non-zero r-derivatives of the 
metric, and these derivatives are the relevant terms in the Einstein equations when 
looking for a solution. So, given the boundary conditions at r = 0, eq.©, we would 
like to obtain the values for the flrst derivatives of the metric at r = e in terms of 
the brane EMT ^ . For doing this we follow the standard procedure of integrating the 
equations of motion in the region with r < e. Notice that if we are to flnd a result that 
is independent of the inner structure of the brane, the dominant terms of the integral 
should be total r-derivatives. If this is the case the value of the integral depends only 
on the value of some functions on the boundary, r = e (and the origin, r = 0), and not 
on the precise solution inside the brane. 

The set of equations we will have to deal with, for the metric at question, will be 
given next. We will offer this equations in a form that makes transparent which terms 
can be integrated exactly and which ones should be neglected when dealing with the 
matching. The iiv components of the Ricci tensor can be written as, 

^LR':; = \[^LK^^' + ^LR^M - v^V^V.L, (4) 

where K^j^y = dj-g^u (we will also use K = Kp, denote four-dimensional {i.e. with 
respect to the metric Qf^y) covariant derivatives and R^y{g) is the Ricci tensor for the 

■^We call brane EMT to TMN{r < e), including possible contributions from the bulk EMT inside 
the extension of the brane. 
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four- dimensional metric g^i^. The 66 component of Ricci tensor reads 

^Li?^ = i^L']' - y^WVpL. (5) 

We see that these two equations can be written as a total r derivative plus terms 
involving only derivatives with respect to the longitudinal coordinates {x'^). As we 
will see in a moment, they will determine the matching conditions. As for the other 
two non-vanishing components of Einstein equations for our metric, the rr and /xr 
components of the Ricci tensor read, respectively. 



and 

R,r = + I^K; + lv^(ir,. - g,.K). (7) 

We can now integrate the fiu and 66 components of the Ricci tensor in the region 
< r < e (the integration in 6 is trivial) to find the desired matching conditions. We 
start with the fiu components, eq.(jl]). Integrating this equation and neglecting the 
terms that do not have r-derivatives'^ one gets 



Jo ^ M, Mt^\Jo Jo ^ 



V9l 



2-K re r 1 1 T'M _ lA/^T 



where we have defined L(x, e) = 1/Mf, (notice that 1/M;, ~ e) and it is understood in 
here and in the following that the subscript |p means that the corresponding function 
is evaluated at r = p. We have also defined the 4D brane EMT, T^, as the integration 
of the full 6D EMT in the region with r < e: 

T^Z^^ dO 1^ dr y^LTj^, (9) 

and T = Tj^J as its trace. 

Now we can repeat the procedure with the 66 equation. Performing the correspond- 
ing integration and neglecting again the terms with only longitudinal derivatives we 

^For the cosmological case we will estimate R^{g) and the other neglected terms of eas. H4l5|l that 
involve /Lt-derivatives, but to obtain the cosmological evolution we need first the matching conditions. 
We will see at the end that these terms are indeed negligible (with respect to the terms we have kept) 
in the cosmological solutions we consider. This is not surprising since once one imposes the constraints 
of having a realistic late-time cosmology the 4D curvature or /x-dependence of the solution have to be 
extremely small. 
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get 

^ r dr J^LRl ^2n(f3- ^] ^ -^-^ dO f dr ^L[T^--T] ^ ~ f , 

VaUJo ^ ' V VgU) Mt^\Jo Jo ^ ' a ^ Mt ' 

(10) 

where we have defined (3{x) = L'{x,e). This equation, together with eq.®, determine 
the exterior space-time geometry associated with a particular energy-momentum stored 
in our codimension two defect. It is apparent that the 66 equation fixes the deficit angle 
in the transverse space, while the /iz/ equations have the clear interpretation of requiring 
a non-zero extrinsic curvature at r = e unless T^ — |5^T = 0. The trace of this quantity 
is referred to as the Tolman mass in 4D cosmic string literature (see e.g. [22]), and it 
is zero for a pure tension brane. For obtaining these equations we have only neglected 
the terms that do not involve r-derivatives in the integrals. 

Notice that the fiu matching conditions are very similar to those of a codimension 
one brane, but because L{x, 0) = we cannot satisfy this equation with a finite 
in the thin brane limit {Mb —>■ oo) in general, in contrast with the codimension one 
case that has a well defined thin limit [18,23]. In fact it is instructive to compare the 
codimension two case with the more familiar codimension one case in more detail, and 
point out their similarities and differences. For a 5D metric like 

ds^ = g^y{x.,r)dx'^dx^ — dr"^, (11) 

we can write the components of the Ricci tensor as 

VaK = \ IVam' + VaKia), (12) 

where, as before, K^u = drQf^u and Ry{g) is the Ricci tensor for the 4D metric gfj,u. We 
can integrate now this equation from r = — e to r = e, a region where we assume that 
some energy-momentum density is localized. We get then 

M^, is now the 5D fundamental mass, is again the integration of the full 5D EMT 
in the (— e, e) region, 

f^^^^Jjr^gn,, (14) 

and T its trace. It is clear in this case that we can take the limit e — keeping T^^, 
and Ry{g) finite if we accept a discontinuity in the extrinsic curvature (notice that 
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we still have y/g\±e \^\o in this limit). We obtain in this way the so called Israel 
matching conditions [23] 

^m! [kX- = - l^'uf, (15) 

since only the extrinsic curvature contribution survives in the thin limit. One could 
think that the analogy between the thick codimension two brane and the codimension 
one case is not surprising since, once the brane has been given a certain width, the wall 
defining the brane is indeed a codimension one object. But there is some information 
in our codimension two matching conditions showing that our system is different from 
a codimension one brane^. First, we have an extra dimensional matching condition, the 
deficit angle contribution eq.()10p. As we have mentioned, in the case of a codimension 
two pure tension brane with = Tg = 0, this contribution absorbs completely the 
effect of the brane on the background, allowing us to keep zero extrinsic curvature. 
This is not the case for the codimension one brane, since the right hand side of eq. ()13p 
does not vanish for a pure tension brane. This difference is what makes codimension 
two braneworlds attractive as a possible solution to the cosmological constant problem. 
Also, in the codimension two case we cannot take the thin limit because VG = Ly/g 
vanishes at r = 0, and then the left hand side of eq.® goes to zero as e — > (even 
allowing for discontinuities in the first derivatives of the metric) if we insist on keeping 

finite. It is then natural to expect that, keeping e finite but small, in a codimension 
two brane-world situation the extrinsic curvature contribution is still the main contri- 
bution to the integral (jH)). Thus, we can safely neglect the integration of the terms 
involving only //-derivatives in most situations, most notably if we are interested in 
late cosmological times for which 4D curvatures are extremely small. 

Another interesting feature of our matching conditions is that we can now use them 
in the fir equation evaluated at r = e to obtain an energy-momentum conservation 
equation for the brane EMT 

The fir component of the bulk EMT determines the fiow of energy- momentum from 
the brane into the bulk. However, even when that term is zero, there can be an 

^This difference comes uitimateiy from tfie boundary conditions we are imposing at r = 0, eqs.©. 
If we were imposing instead tlie boundary conditions L'{x,0) — and L(x,0) different from zero, we 
would be describing a 4-brane with a compact dimension (9) and a Z2 symmetry at r = 0. 
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exchange of energy-momentum between the longitudinal and the transverse directions 
(on the brane) correlated with a possible space-time dependence of the brane width 
that could be interpreted by 4D observers as apparent violations of the conservation 
of energy-momentum. 

Up to now we have performed an integration of the 6D Einstein equations in the 
space-time region filled by the brane in order to obtain the matching conditions. We 
have identified a set of terms in the equations that can be integrated in a model 
independent way, and we have seen that these terms are indeed the dominant ones 
provided the /x-dependence of the solution is small. In this way we have obtained 
eqs. ()8ll0p that relate the first (transverse) derivatives of the metric just outside the 
brane with the integrated brane EMT. We are taking g^v\t as our "induced metric" 
for the defect, since we are evaluating most functions in r = e when dealing with 
the matching. In fact we see that all the functions appearing in the matching are 
evaluated at the brane boundary except for the ratio ^/g\o/ ^/g\e appearing in the 00 
matching condition, eq. ()10|) and in the EMT conservation equation, eq. ()16p . In order to 
avoid making reference to functions evaluated inside the brane when dealing with the 
matching we will consider that this ratio can be approximated by one in this equation. 
For this we just need that 



« 



rrifi 



(17) 



This is because we can put a bound on the difference y/g\^ — ^/g\o as 

V9\e - V9\o < v^lo — Tr~ ~ v^lc 



2Mb 

where K\max is the maximum value of the function for r < e and we have used = 
y/gK/2. Using now the matching conditions, eqs. ()8pi()|l . we arrive at the mentioned 
requirement, eq. (fT7|) . When this condition holds, we could speak of a quasi-pure- 
tension-brane, and the extrinsic curvature (times the brane width) is negligible with 
respect to the deficit angle. This allows us to approximate the ratio of the metric 
determinants in eq. ()l()|l by one. We quote here the actual matching conditions again, 
assuming such condition holds 

27r(l -(3) - -}-(-f- fi] ~ -^rt^ (19) 

(20) 



M, 



b 
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We would like to emphasize the generality of these matching conditions. They apply 
to any codimension two brane, provided we can neglect the longitudinal /i-derivatives 
when compared with the transverse r-derivatives in the solution inside the brane, 
and our condition eq. lfTTj) holds. For instance, it is now straightforward to interpret 
several known solutions with naked codimension two singularities as being sourced by 
a codimension two object with certain energy-momentum. Consider as an example 
metrics that near r = can be approximated as 

Qfiu ^ K,ir°''ViMu + • • • , (21) 
L ~ K2r''^ + .... (22) 

Our matching conditions then imply 

J^2a^ ^ _To_ . . 

where we have taken = Tq 5^ and we are assuming that eq. (jl7|) holds. We see how 
the required energy-momentum for the defect sourcing these solutions depends on the 
brane width. The brane width acts as a cut-off for the curvature, and gets rid of the 
singularity once one considers a thick defect: remember that these relations have been 
obtained matching a regular geometry at r = (that implies eqs.Q) with the exterior 
geometry given by eqs. (j21l22j) . It is interesting to point out that in the thin brane 
limit {Mh oo), the brane tension diverges if 02 < 1 while it goes to zero if 02 > 1 
(in this latter case we can not satisfy our assumption in the thin limit, eq. (jl7p . and 
the first matching condition above would have some corrections). We can recognize 
the only parameters with a well defined thin brane limit yielding a finite brane energy- 
momentum (ai = 0, 0^2 = 1) as a purely conical geometry. In particular the singular 
solutions of 6D supergravity found in [13], can be cast in the form given by ecs. ()21l22|) 
and interpreted as being sourced by a codimension 2 defect with an energy momentum 
tensor given by the formulae above. 

Also, we can now match an exterior AdS geometry with a regular geometry on r = 
and check what kind of energy-momentum distribution supports such spacetimes. This 
kind of exterior geometry has been obtained as the spacetime produced by the Nielsen- 
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Olesen vortex of the Abehan Higgs model in 6D^ [25] which, interestingly enough, 
localizes not only gravity but also gauge interactions at the vortex core [26] (see also [27] 
for the non-abelian case). Applying our matching conditions to an exterior geometry 
like (AdSe space with a compact dimension) 



ds' = e^'''r]^,dx''dx'' - Lle^'^'de'' - dr\ (25) 



where k = ^ — ^i, and A is the bulk cosmological constant, we obtain the required 
brane EMT as 

Lok ±k To 

e ~ — — , (26) 



±^ = -ArAT^ + Tn, (27) 



2ttM^ ' 
1 

"Mb ~4^4 

where we have taken again = Tq 5^ and we are assuming that eq. (fT7j) holds. Notice 
that if we choose the minus sign for ±k, so the volume of the spacetime is finite, 
we need Tq > 27rAf^. This condition does not mean that the curvatures are in the 
solution bigger than the fundamental mass, since we expect that ~ ^/T^^ ~ T^/^h- 
For low values of Mf, we can still have the 4D brane EMT (T^) of the order of 
or bigger while the 6D one (T^) is hierarchically smaller and the curvatures in the 
full 6D solution are under control, but one can not keep curvatures under control in 
the thin brane limit (see also the discussion in [25]). However, if we want to restrict 
ourselves to weakly gravitating branes we should choose the plus sign, and then the 
extra dimensional volume is infinite. This is because L is a growing function at the 
origin with a positive r-derivative, and then it is "easier" to match the geometry with 
an exterior L function that also has a positive derivative (and this chooses the plus 
sign above). 

Our matching conditions relate brane parameters with metric deformations in its 

surroundings, but they do not tell us much about the phenomenological viability of 

these braneworlds. In an ideal situation, for a given brane energy-momentum and 

width, and for a given bulk EMT we would impose the matching conditions in the bulk 

solution (as a perturbation, perhaps, over a known static solution) and find out the 

implications for the brane induced metric, g^y{x^e). The relation between the brane 

^For the global vortex the exterior geometry is AdS^ x [24], but this case is not a solution for 
a pure cosmological constant in the bulk. 
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induced metric and its EMT would determine what type of gravity brane observers 

would feel. But this approach is usually very difficult to implement in practice, since 

it is very hard to find analytic bulk solutions (even perturbative ones, see [16] for work 

in this direction). It is however possible to obtain a good deal of information on the 

curvature of the brane induced metric without actually solving the bulk equations. The 

idea, that we will elaborate on in more detail in the next section for the cosmological 

case, is to evaluate Einstein equations just outside the brane. In particular, the rr 

component of the Einstein tensor, evaluated at r = e reads, 

1 A/f2 
G;|. = -2fi(9) + ^P'; + T»] 

^ '-[4T^T^ - {TPf - 2t;(t; + T's) - (t; + r,^)(5T; - 3t,^)] - ^ 



*(28) 

where we have used the matching conditions eqs. p9l20|) and have neglected the term 
MfoVVpLle. It is not clear a priori that this term is negligible as compared with 
the induced metric curvature (the integrals of both terms have been neglected in the 
matching conditions). However it is in general related to /i— derivatives of (3 that, 
through our matching conditions, can be argued to be negligible under the assumptions 
we are using. Again, we will be more specific about the size of the different terms 
when discussing cosmological solutions. This equation determines the curvature for the 
induced metric in terms of the brane EMT and the bulk EMT evaluated at the brane 
boundary. This equation will become much more illuminating when particularised to 
the cosmological case as we will see in the next section. 

3 Cosmology on a codimension 2 brane. 

In order to study the cosmological implications of our matching conditions we partic- 
ularise the metric ansatz to 

ds'^ = N{t, rfdt^ - A{t, rfdx^ - dr^ - L{t, rfdO'^, (29) 

where we have taken flat spatial sections in the brane for simplicity and we set A^(t, 0) = 
1 by performing a redeflnition of the t coordinate. In this section we will assume 
that eq. (jl7j) holds in our system, so we can use the simpler version of the matching 
conditions, eqs. (jl9l2(J|) . instead of the more general one, eqs. (j8llUp . Remember that. 
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since we are assuming that the first derivatives are small compared with the brane width 
we also have, at the level of approximation we are working, N{t, e) ~ N{t, 0) = 1, and 
also ^4(^,0) ~ A{t,e) = a{t). So the matching conditions, eqs. ()19l20|) . take now the 
form 

A'L Mb . ^ . . 

N% = ^^MP + P)+Pr+Pel (31) 
^^-f^^ = -^P-Pr + ^Pe]- (32) 

where we have taken a "cosmological" brane EMT: T^^ = diag{p, —p, —p, —p, —pr, —pe)- 
These equations, however, do not yield any information about the cosmology one can 
expect in these models. One constraint on it is given by the brane EMT conservation, 
eq.(fTB|). that reads for our cosmological set up, 

P + 3- P + P - dtPr - [Pr +Pe)^ = 27r^!^. 33 
a Ml) Mb 

However, we would like to obtain the equations that govern the evolution of the scale 
factor of the brane, a{t), in terms of the brane EMT. As we said at the end of the 
previous section, with our matching conditions we could, in principle, find the bulk 
solution for a given model (as a perturbation, perhaps, over a known static solution) 
and figure out the implications of this perturbed solution for the time dependence of 
the scale factor. Instead of doing that we will show that a lot of information on the 
cosmology of codimension two branes can be obtained by evaluating Einstein equations 
just outside the brane, following the same approach as [11] for the codimension one 
case. The matching conditions tell us what the first r— derivatives of the metric are 
at the brane boundary, eqs. (|3(JH3^ . whereas second r— derivatives can only be found 
by solving Einstein equations. The crucial point to note is that, out of the five non- 
zero (for the cosmological metric) components of the Einstein tensor, the tt, xx and 66 
involve second r-derivatives and therefore allow us to algebraically compute their values 
at the brane boundary whereas the tr and the rr components do not involve second 
r— derivatives but only known first r— derivatives and first and second time derivatives 
of the metric. The former gives the brane EMT conservation, eq.(|33p whereas the 
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latter, evaluated at r = e, reads 




a 

+ - 
a 



[3(p + pf + {pr + pe) [2(p - 3p) - + 3pe 




(34) 



where we have particularised eq. ()28|l to the cosmological case. Here we can be a bit 
more specific about the size of the term we are neglecting, which is (we do not explicitly 
write factors of order one) 



where in the third equality we have approximated (as a conservative order of magni- 
tude) dfL'{r < e) ^ dfL'\^, and in the fourth one we have used the matching conditions. 
We therefore see that, at least at late cosmological times it is utterly negligible as com- 
pared with the terms we are keeping. Nevertheless it should be noted that this is just 
an order of magnitude estimation and one should carefully check this approximation 
when dealing with particular models (for which this term could play an important role 
in the cosmology of our brane). The equation we have obtained is a generalised Fried- 
mann equation that incorporates the matching conditions for our general codimension 
two brane. Taken together with the energy- momentum conservation equation, eq. ()33|) . 
suffices to determine the evolution of the scale factor and therefore the cosmology. As 
a first check to this equation we can note that we recover the expected behaviour in 
the case of a pure tension brane with no rr and 66 components of the EMT. This 
case corresponds to p + p = Pr = Pe = and we get that the expansion of the brane 
depends solely on the bulk EMT and not on the brane tension, agreeing with the solu- 
tions presented in the literature for this case [12]. As a matter of fact, we could have 
guessed the generic form of the modified Friedmann equation based on the well defined 
infinitesimal limit for a pure tension brane as something like 



where J-',Q are arbitrary functions with the only restriction that Q{0,0, p,p) = and 
JF(0) = 0. The bulk EMT might also have some implicit dependence on p,p and 
Pr,e- A critical point when considering self-tuning issues is the particular form of the 




(35) 




(36) 
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function T. If it was linear in p + p it could represent an important step towards a self- 
tuning scenario yielding a cosmology of the type studied in [28] whereas the quadratic 
dependence we have actually found would be in conflict with phenomenology. But this 
negative conclusion in the self-tuning issue is of course a bit premature, and should not 
be taken too seriously, in the sense that we have not included the possible dependence 
of the bulk EMT on the brane parameters. A reahstic model might exist that has 
self-tuning features hidden in such dependence. For the time being however we will 
naively assume that the bulk EMT does not have any implicit p, p or pr,6» dependence, 
and we have just a cosmological constant in the bulk. In that case, it is possible to 
obtain a realistic cosmology if one is ready to give up self-tuning considerations. In our 
universe, the term in this equation proportional to (p would be extremely small 
(~ Pmatter)^ whilc (assumiug pr and pe are constant) the terms that go like pr + pe or 
iPr + Pe){5pr — 3pd) would act as a cosmological constant. One expects then that the 
term cx {pr + Pb){p — 3p) would give the dominant time-dependence and therefore a 
conventional cosmology 

where we have just neglected the terms proportional to (p + and p^ Q and defined 
T^\^ = Ar. There remains, of course, the issue of the effective cosmological constant 
that would have to be tuned to zero. We have several parameters that can be chosen 
at will in the equation above, so one can fine-tune the effective cosmological constant 
to a small value by requiring 

where we have considered that p ~ —p ~ Tq. There is a priori no reason to expect 
such a cancellation, so this scenario does not seem to yield any light on the cosmolog- 
ical constant problem. But notice that the value of pg and Pr coming just from the 
integration of the bulk cosmological constant inside the brane has an order of mag- 
nitude that, for a weakly gravitating brane (Tq/M^ << 1), already agrees with this 
requirement: in case of having no "brane" contributions to these quantities one would 
expect Pr,Pe ~ — 27rA,.e^ ~ —27rAr/M^. 

It is very interesting that such a simple model can yield a realistic cosmology, but 
one would expect that in more involved models the bulk deformation produced by the 
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brane EMT would affect also. Depending on how reacts under a deformed so- 
lution (satisfying our p- and p-dependent matching conditions), one obtains different 
cosmologies. We see that considering a constant = Ar, Pr and pe can yield a con- 
ventional cosmology (with the cosmological constant problem included), but it would 
be very interesting to see what happens in other cases. In the case of having two com- 
pact extra dimensions stabilised with a magnetic flux for instance, the deficit angle in 
the transverse space affects the local energy density of the flux, and therefore T^\e in 
eq. ()34|) . and this could be the dominant effect when determining the cosmology [16]. 
A deeper study of these issues is currently under way. Nonetheless, it is interesting 
to point out that in the case of having just a cosmological constant in the bulk, the 
effective Planck mass is 

Mh = -8-^^(P. + Per - 4vr^ fl - , (39) 



where we have used eq. fl38j) . Any relation with the extra-dimensional volume is not 
transparent in this equation but let's imagine that we are in a situation whose exterior 
metric is approximately described by the geometry given by eq.(|25|). simply AdS space 
with a compact dimension, and the matching conditions reduce to eqs. (j26l27|) . Taking 
the minus sign in (j25|) . so the extra dimensional volume is finite, we can write 

V2 ~ 27r r Loe'^'^'/'dr = (40) 
Jo 5k 



We can use now the matching conditions and the relation k = y — |^ to get 

Mh = -jMtV,. (41) 

We need to chose the negative sign in the metric ()25j) in order to get a positive Planck 
mass squared. If we want to match our brane with an exterior AdS geometry as in 
[25] , and obtain a realistic cosmology, we need that To > 27rM^ (that do not necessarily 
imply curvatures of order ~ M^, as we have previously commented). If on the other 
hand we insist on a weakly gravitating brane, Tq << M^, then our matching conditions 
show that the appropriate branch is the plus sign in (|25|) and so the volume of the extra 
dimensions is infinite. Interestingly enough, the effective 4D Planck mass remains 
finite although imaginary, a situation identical to the codimension one case with a 
negative tension brane [10]. Also, as we will see below, the dependence on the different 
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scales with high exponents makes it very easy to generate a hierarchically large Planck 
mass without departing from our approximations. One can point out again several 
similarities between this possibility and the well known codimension one braneworlds 
consisting on a 3-brane moving in a 5D AdS space [10,11]. In both cases the brane 
EMT enters in the definition of the effective Planck mass, and in both cases we need 
a non-zero "brane vacuum energy" (the brane tension in the codimension one case, 
or the Pr + pe parameter in our case) in order to find a realistic cosmology. This 
is a result of the quadratic dependence of the generalized Friedman equation on the 
brane EMT parameters in both cases. Also, in both cases we need to fine tune the 
bulk cosmological constant against brane parameters in order to obtain a small effective 
cosmological constant on the brane. It is also thanks to this fine-tuning that we recover 
the relation of the effective 4D Planck mass with the higher dimensional one times the 
extra-dimensional volume, as expected from KK arguments (as in the codimension one 
case). It would be interesting to push this analogy further, and we might be able to 
interpret our thick codimension two brane as a "curled up" codimension one brane, 
since as we have commented previously, the wall defining the brane boundary is indeed 
a codimension one hypersurface. 

We can also check now the magnitude of the terms we are neglecting in the Einstein 
equations when doing the matching for these solutions. In order to do so, we consider as 
an example particular values of the different parameters, motivated by a "TeV brane" . 
We consider a weakly gravitating brane situation, and in this case the exterior geometry 
would not be of the type given by eq.(j2Sl), since we need a positive bulk cosmological 
constant in order to obtain a realistic cosmology (see eq. ()39|) ). Taking for instance the 
brane parameters to be of the order 





Mb^ TeV, 



(42) 
(43) 
(44) 



and the fundamental scale 



M, ^ 1.7 X 10^ TeV, 



(45) 



so that the Planck mass comes in the right size 



Mpi 



lO^^GeV, 
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we see that the terms we neglected in the matching are of order 

/' * ^ ~ m4^^' + 3P - Pr + 3p„) ~ 10-'». (47) 

We have assumed that the rr component of the bulk EMT is fine-tuned to give a small 
(realistic) Hubble parameter, aI^^ ~ 200 GeV and for the second equation we have 
used arguments similar to the ones leading to eq. It is therefore clear that these 

terms are indeed smaller than the ones we have considered, that are of order 

^ drL" = 1 - ~ _Lj(p - 3p + p, - 3pe) ~ 10'\ (48) 

r driLK)' = # ~ - 10--. (49) 

Note that these terms, although much larger than the neglected ones, still have a 
hierarchy between themselves, as prescribed by our requirement, eq. p7|) . As we said 
we need the small cosmological constant fine-tuning so that the contributions given by 
fl46j) are smaller than the ones coming from the terms like P9|l . One can understand the 
hierarchy of the different terms noticing that the terms we neglected in the matching 
are in these solutions proportional to the brane EMT squared or its time derivatives 
(assuming that the term linear in + pq in eq. ()37|) is cancelled by A^), while the 
terms we have kept are matched with the brane EMT linearly. One should, however, 
check for any particular solution the level of approximation that using our matching 
conditions represent, since it is not guaranteed in general that they constitute a good 
approximation. In case they do not, our matching conditions have corrections that 
depend on the internal structure of the brane. A numerical estimation of these terms 
would be advisable when using our matching conditions, although as we have seen one 
can easily find models with solutions in which they are utterly negligible, and therefore 
in these models the use of our matching conditions is fully justified. 

Before finishing this section we would like to stress that this particular example (and 
its associated cosmological constant fine-tuning) has nothing to do with a possible self- 
tuning mechanism. In these models, solutions in which the brane is curved exist on 
equal footing with solutions in which the brane is fiat. So the required tuning should 
be provided by extra considerations, the most promising being supersymmetry (see 
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[29] and references there in). The crucial point is that a thick brane, via our matching 

conditions, allows for a dynamical deficit angle and therefore the possibility of self- 
tuning is present in such a set-up. A more detailed study of self-tuning issues in these 
models requires knowledge of the bulk solution in more involved models and is currently 
under way. 

4 Conclusions, open issues and future prospects 

The goal of this letter was to study the dynamics of the induced metric on a codimension 

two braneworld, and in particular we were interested in its relation with the brane 

EMT. The solution is singular in general for an infinitesimally thin brane^, except for 

the case of a pure tension brane. Contrary to what might seem, this case is not relevant 

for the study of the so called self-tuning properties on codimension two braneworlds 

due to the staticity (for the deficit angle) of the solution as opposed to the intrinsically 

dynamical nature of self-tuning. Due to these reasons we had to abandon the thin 

brane idealization and consider a brane of finite thickness. Our first step was to find 

the matching conditions in the second section, i.e., the set of equations that relate the 

brane EMT with the deformation of the surrounding spacetime it produces. Since we 

are dealing with a thick brane that has a singular thin limit, it was not clear that one 

should be able to do this without knowing all the precise small-scale structure of the 

brane. We have however shown that one can obtain this set of equations depending only 

on the integral of the brane EMT as a good approximation provided certain conditions 

are met. These conditions have the interpretation of requiring that the brane has an 

energy-momentum lying mainly along the parallel directions (so it is close to a pure 

tension brane), with small 4D curvature. We have also seen in this second section that 

the fir component of the Einstein equations gives, when evaluated just outside the 

brane using our matching conditions, the energy-momentum conservation equation for 

the brane. Using our matching conditions we have been able to interpret some singular 

solutions of 6D supergravity found recently in [13] as being sourced by a codimension 

two defect with certain energy-momentum and we have also paid particular attention 

we have previously mentioned, this is not the case when one considers Einstein-Gauss-Bonnet 
gravity on the bulk [15], and in fact one can find non-singular solutions even for infinitesimally thin 
higher codimension braneworlds when one uses the general Lanczos-Lovelock Lagrangians in higher 
dimensions [21]. 
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to the simple case of having just a cosmological constant in the bulk. In this case 
wc have matched our brane with an exterior AdS geometry as in [25], obtaining the 
required brane properties. 

In the third section, attempting to obtain further information that allows us to 
assess the phenomenological viability of these models, we have specialized our equations 
to the cosmological case. The matching conditions relate however the brane EMT with 
the first derivatives of the metric with respect to r, the orthogonal coordinate, and to 
obtain the cosmology we would hke to relate the brane EMT with the parallel (time) 
derivatives of the metric. Fortunately, this can be done simply by using the matching 
conditions in the rr component of the Einstein equation evaluated just outside the 
brane, at r = e. The reason is that this component of the full 6D Einstein equations 
is the only one (apart from the /xr) that does not involve second r— derivatives, so we 
can get rid of the first ones using our matching conditions and we are left with only 
time derivatives of the induced metric and brane parameters, obtaining the modified 
Friedmann equation we were looking for. This procedure is completely analogous to 
the one followed in the codimension one case by Binetruy et al. in [11], the only 
difference being the added complication here of having to consider a thick brane, since 
the thin limit is singular for the codimension two case. We have identified a model that 
could yield a realistic cosmology, simply considering a constant value for the orthogonal 
components of the brane EMT and fine-tuning the bulk cosmological constant to get 
a small effective cosmological constant for the brane metric. Put like that, this model 
docs not seem to shed any light on the possible self-tuning behaviour of codimension 
two brane worlds, one of the original motivations for considering these class of models. 
It would be interesting in any case to explore this simple model with just a cosmological 
constant in the bulk in more detail, in particular to study which bulk geometries are 
obtained from it when we have a positive cosmological constant in the bulk, that as we 
saw, makes compatible a weakly gravitating brane with a realistic cosmology. In case 
of having a negative bulk cosmological constant, the exterior geometry is in some cases 
just AdS space with a compact spatial dimension [25] , and for a pure tension brane the 
bulk is simply a wick rotation and analytical continuation of the AdS-Schwarzschild 
geometry (see below), so one might expect that the solution asymptotes to AdS space 
in general when we have a negative cosmological constant in the bulk. 

But in fact the main uncertainty, and source of model dependence in our cosmolog- 
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ical equations is the bulk EMT: our modified Priedmann equation depends on its rr 
component. One can expect that in particular models the T/vfAr— dependent matching 
conditions for the brane imply that this bulk EMT component also gets an implicit Tmn 
dependence, i.e. T^\e would also depend on p, p, p,. and pe for the cosmological case. It 
is certainly conceivable that this dependence is the dominant one, and one could then 
obtain different cosmologies depending on the particular model, or compactification, 
one is dealing with. A more careful examination of these issues for different models 
will be deferred to a future publication. It is worth pointing out that the bulk and the 
brane curvatures are independent parameters that in principle have no reason to be 
related, even for some given brane parameters. This can be seen by considering the 6D 
black hole with cosmological constant [30], substituting in it t — >^ i9 and analytically 
continuing the constant (r, 9) hypersurfaces to a Lorentzian manifold of curvature H^. 
Then one can interpret the solution as a pure tension infinitcsimally thin codimension 
two braneworld where is the curvature of the brane induced metric [31]. One can 
see then that the brane curvature, the bulk cosmological constant (that determines 
curvature of the bulk) and the brane tension are independent parameters of the solu- 
tion. So these models do not have any selftuning behaviour per se, if certain models 
only admit flat 4D geometries [13], supersymmetry is to blame, and supergravities with 
these properties are also known in 7D [32] (where there are 3 compact dimensions and 
4 fiat dimensions). The problem, of course, is to obtain an effective 4D theory with 
supersymmetry broken at a high scale without spoiling 4D flatness, and codimension 
two branes could help in this [29] . It would be very interesting to impose our matching 
conditions in the solutions of the 6D supergravity that have been proposed in order 
to realize the selftuning behaviour [13], and in particular to study the implications of 
these matching conditions for the supersymmetry of the background, but as we said 
before, a closer examination of these issues is beyond the scope of the present paper, 
and will be deferred for future work. 
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